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COMPARISON  WITH  EXPERIMENT  OF  SEVERAL  METHODS 
OF  PREDICTING  THE  LIFT  OF  WINGS  IN 
SUBSONIC  COMPRESSIBLE  FLOW 
By  Harry  E.  Murray 

SUMMARY 


Several  methods  of  predicting  the  lift  of  wings  in  subsonic  com¬ 
pressible  flow  were  compared  with  experiment.  An  experimental  verifica¬ 
tion  of  Kaplan* s  formula  for  the  effect  of  compressibility* on  the  lift 
of  wing  sections  was  obtained. 

Semiempirical  formulas  were  developed  for  predicting  the  subsonic 
effects  of  compressibility  on  the  lift  of  finite-span  wings  based  on 
corrections  to  the  section  lift-curve  slope.  These  semiempirical 
formulas  yielded  better  agreement  with  experiment  than  previously 
derived  theoretical  methods.  The  agreement  at  small  sweep  angles  was 
slightly  better  when  thickness  was  considered  in  the  semiempirical 
formulas . 

.  Both  experiments  and  calculations  indicated  a  decrease  in  the 
variation  of  lift  with  Mach  number  for  increasing  sweep. 


INTRODUCTION 


The  effect  of  compressibility  on  the  lift  of  finite— span  wings  has 
been  extensively  discussed  in  previous  papers  (references  1,  2,  and  3, 
for  example).  These  papers  discuss  compressibility  effects  in  terms 
of  an  affine  transformation  based  on  small— perturbation  theory  herein 
referred  to  as  the  three-dimensional  Prandtl  transformation.  The 
application  of  the  three-dimensional  Prandtl  transformation  to  the 
lifting-line  theory  of  unswept  wings  is  discussed  in  references  1  and  2. 
An  application  of  Weissinger*s  approximate  lifting-surface  theory  of 
wings  of  arbitrary  sweep  (reference  4)  is  discussed  in  reference  3. 

When  compared  with  experiment ,  these  existing  methods  did  not  yield 
entirely  satisfactory  results. 

Kaplan  (reference  5)  has  shown  that  including  the  thickness  of  a 
two-dimensional  airfoil  in  calculations  of  the  effect  of  compressibility 
on  the  lift  results  in  appreciable  effect  at  high  subsonic  Mach  numbers. 
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Inclusion  of  thickness  in  the  three-dimensional  case  consequently  may 
also  have  an  appreciable  effect.  Eather  formidable  mathematical 
difficulties  are  encountered  in  any  rigorous  attempt  to  consider  the 
thickness  of  a  finite— span  wing  in  subsonic  compressible  flow.  To 
attempt  an  approximate  adaptation  of  Kaplan* s  two-dimensional  solution 
to  finite— span  wings  therefore  seems  reasonable.  In  the  present  paper 
such  an  approximate  adaptation  of  Kaplan*s  results  to  swept  and  unswept 
wings  is  obtained.  The  method  is  based  on  lifting-surface  theory  (refer¬ 
ence  6).  Available  test  data  are  compared  with  the  present  method,  as 
well  as  with  several  other  methods  which  neglect  the  effect  of  airfoil 
thickness. 


SYMBOIS 


finite— span-^wing  lift  coefficient  — 


£cv2 


CL 

a 


5a 

5a 


M 


A 


section  lift  coefficient  /^ift  span 

'  2 

air  mass  density 
wing  area 
wing  chord 
free^tream  velocity 
angle  of  attack: 


section  lift-curve  slope  at  M  =  0  1 

section  lift-curve  slope  at  Mach  number  of  M 

/  Y  \ 

Mach  number  ( - j 

\Velocity  of  sound/ 

sweep  angle  of  wing  quarten-chord  line,  positive  for  sweepback 


Me  =  M  cos  A 


M-  = 


/l  -  M2 
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Kaplan  (reference  5)  has  presented  a  method  whereby  the  effect  of 
airfoil  thickness  can  be  considered  in  the  calculation  of  the  lift  of  a 
two-dimensional  wing  at  an  angle  of  attack.  Kaplan* s  formula,  for  the 
effect  of  compressibility  on  lift  is 

^  =  IT=^  +  |(l-  e~2^  ^  “  1)  +  +  1)(^2  ~  l)2]  (1) 

The  effect  of  thickness  is  included  in  this  formula  by  means  of  the 
parameter  \  which  relates  the  chord  of  the  airfoil  to  the  radius  of 
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the  conformal  circle  in  Theodor sen’s  potential  theory  of  arbitrary 
profiles  (reference  7)*  As  the  airfoil  thickness  approaches  zero,  \ 
approaches  zero  and  Kaplan’s  formula  can  he  seen  to  approach  the  Prandtl 
(or  G-lauert)  factor,  which  is 


(2) 


The  variations  of  section  lift-curve  slope  with  Mach  number  for  the 
NACA  66,1— 115  airfoil  and  for  the  same  airfoil  with  a  beveled  trailing 
edge  (models  1  and  2,  table  i)  are  shown  in  figures  1  and  2,  respectively 
(from  reference  8).  Calculations  applying  Kaplan’s  formula  (equation  (l)) 
and  the  Prandtl  factor  (equation  (2))  are  also  shown  for  comparison. 

The  short  vertical  lines  on  the  curves  of  figures  1  and  2  indicate 
the  lower  limit  of  the  test  data  below  which  the  experimental  curve 
is  extrapolated  to  zero  Mach  number.  Consideration  of  the  thickness 
(the  Kaplan  method)  improves  the  agreement  between  theory  and  experiment . 
Similar  agreement  is  shown  in  figure  3  for  'an  airfoil  approximating  the 
NACA  0012—64  airfoil,  designated  as  E— 0009  in  reference  9,  at  zero  sweep 
(velocity  of  free  stream  normal  to  wing  span) . 

In  reference  10,  Jones  indicated  that  the  effect  of  compressibility 
on  a  swept  two-dimensional  wing  is  the  same  as  on  an  unswept  wing  in  a 
stream  of  reduced  Mach  number  Me,  where 

Mg  =  M  cos  A  (3) 

Kaplan’s  formula  can  thus  be  adapted  directly  to  the  two-dimensional 
swept  wing  by  replacing  u  by 


to 


1/1  -  %  /1  _ 


(4) 


M2cos^A 


Calculations  made  with  this  modification  of  Kaplan’s  formula,  together 
with  experimental  data  and  calculations  made  with  the  Prandtl 
correction,  are  presented  in  figure  3  for  "the  two-dimensional 
NACA  0012-64  airfoil  (model  3,  table  i)  at  three  angles  of  sweep. 

The  results  indicate  that  the  foregoing  method  of  predicting  the  effect 
of  compressibility  on  a  two-dimensional  swept  wing  produces  good  agree¬ 
ment  with  experiment  for  A  =  20°  when  thickness  is  considered.  The 
reason  for  the  somewhat  poorer  agreement  for  A  =  40°  is  not  known. 

Below  the  short  vertical  lines  on  the  curves  of  figures  1,  2,  and  3 
the  experimental  curve  is  extrapolated  to  zero  Mach  number  by  means  of 
Kaplan’s  formula  in  which  is  used. 
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The  Lift  of  Three-Dimensional  Wings 

The  effect  of  compressibility  on  the  lift  of  three-dimensional  wings 
is-  calculated  by  four  different  methods.  All  four  methods  stem  from  the 
three-dimensional  Prandtl  transformation.  According  to  the  three-dimensional 
Prandtl  transformation  as  set  forth  in  reference  3,  the  effect  of  com¬ 
pressibility  on  lift  can  be  obtained  by  increasing  by  the  factor  —  - 1 

\fl  -  M2 

the  lift  for  incompressible  flow  of  wings  having  equivalent  aspect  ratios 
given  by 


Ac  =  A/l  -  M2 


and  equivalent  sweep  angles  given  by 


.  .  tan  A 

tan  A  = 

yr  -  m2 


Methods  1  and  2  are  based  on  an  interpretation  of  the  three-dimensional 
Prandtl  transformation,  which  is  strictly  correct  only  for  unswept  wings 
of  high  aspect  ratio  to  which  lifting-line  theory  is  applicable.  This 
particular  interpretation  was  adapted  for  application  to  wings  of  moderate 
aspect  ratio  with  sweep  because  it  afforded  a  simple,  logical  means 
wherein  a  correction  for  thickness  based  on  Kaplan’s  results  (reference  5) 
could  be  applied.  Methods  3  and  4,  however,  are  strict  applications  of 
the  three-dimensional  Prandtl  transformation.  No  logical  or  practical 
method  could  be  discovered  for  adapting  a  strict  application  of  the  three- 
dimensional  Prandtl  transformation  to  the  purpose  of  accounting  for  the 
effects  of  wing  thickness. 

Method  1.—  According  to  reference  1,  if  the  three-dimensional 
Prandtl  transformation  is  applied  to  lifting-line  theory,  the  lift-curve 
slope  of  unswept  wings  is 


¥ 


where  Ac  =  £A  and 


A^i 

Ac  +  aq  21^2. 


Equation  (5) 


(5) 


also  can  be  written  as 


•4 


A  + 


H  37.3 

(3  it 


(6) 


6 


MCA  TN  No,  1739 


Equation  (6)  shows,  as  has  already  "been  observed  in  reference  2,  that 

for  lifting— line  theory  the  effect  of  compressibility  can  be  accounted 

for  by  simply  using  the  correct  variation  of  section  lift-curve  slope  * 

with  Mach  number.  Equations  (5)  and  (6)  apply  rigorously  only  to  the 

limiting  case  of  unswept  wings  having  very  high  aspect  ratio.  According 

to  reference  6,  the  lift-curve  slope  of  unswept  wings  in  incompressible 

flow  can  be  obtained  more  correctly  than  by  lifting-line  theory  from  the 

following  equation 


A ai 

°La = « 777733 

e  1  it 

which  is  based  on  lifting-line  theory  but  corrected  according  to  lifting- 
surface  theory.  The  product  AEe  is  shown  in  figure  If  the  effect 

of  compressibility  can  be  assumed  to  be  accounted  for  by  correcting  the 
section  lift-curve  slope,  as  in  the  case  of  lifting— line  theory,  then 
equation  (7)  can  be  written  as 


AEe  + 


&i  57^3 
P  * 


(8) 


Equation  (8)  is  not  £  strict  application  of  the  three-dimensional  Prandtl 
transformation  because  of  the  presence  of  the  quantity  Ee.  Tne  effect 
of  thickness  can  be  approximately  introduced  into  equation  (8)  as  a 
correction  to  the  section  lift-curve  slope  by  substituting  p^-  for  P; 
the  resulting  equation  is 


to 


fl  57^3 


(9) 


Reference  11  has  shown  that  the  effect  of  sweep  can  be  approximately 
accounted  for  in  formulas  for  the  lift— curve  slope  of  unswept  finite- 
span  wings  by  multiplying  the  section  lift-curve  slope  by  the 
factor  cos  A.  If  some  new  factor  is  assumed  to  account  exactly 

for  sweep,  then  equation  (9)  can  be  written  to  include  sweep,  as  follows: 


A 


A 


AEe  + 


57-3 

Jt 


i 


(10) 
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For  a  swept  wing  the  results  of  the  two-dimensional  analysis  will  he 
used  to  account  approximately  for  the  variation  of  compressibility  effect 
with  sweep.  For  this  purpose  (Bg-  in  equation  (10)  will  be  based  on  jiA 

or,  in  other  words,  upon  the  component  of  Mach  number  normal  to  the 
swept  panel.  Comparisons  of  equation  (10)  with  experiment  indicated 
that  best  agreement  resulted  when  A  is  defined  as  the  sweep  of  the 
quarter— chord  line  and  X,  as  used  in  the  equation  for  (Bg-,  corresponds 

to  the  airfoil  section  normal  to  the  quarter-chord  line. 

Because  FA  cannot  be  accurately  evaluated  from  available  theories 
for  all  useful  values  of  aspect  ratio  and  sweep,  equation  (10)  is- 
valuable  because  of  its  ability  to  predict  the  effect  of  compressibility 
on  a  wing  for  which  low-speed  experimental  data  are  available.  The 
quantity  a^FA  can  be  evaluated  from  the  low-speed  test  condition  and 

is  assumed  to  be  independent  of  Mach  number  in  the  subsonic  range. 

Equation  (10)  has  been  applied  as  herein  explained  and  is  compared 
with  experiment  in  figures  5  to  19  as  method  1.  The  short  vertical  lines 
on  the  curves  of  figures  5  to  19  indicate  the  lower  limit  of  the  test 
data  below  which  the  experimental  curve  is  extrapolated  to  zero  Mach 
number  by  means  of  method  1. 

Method  2.-  Method  2,  which  is  identical  to  method  1  except  that 
X  =  0,  is  included  in  figures  5  to  19  in  order  to  indicate  the  order 
of  magnitude  of  the  effect  introduced  by  thickness.  When  X  =  0,  pE 
is  replaced  by  |B  in  equation  (10),  and  the  resulting  equation  is 


AEe  + 


aq 

P 


ZLl 

7t 


Method  3.-  A  strict  application  of  the  three-dimensional  Prandtl 
transformation  (reference  3)  to  very  thin  unswept,  or  very  slightly 
swept,  wings  transforms  equation  (7)  to 


Acai 


AcEe, 


+  a, 


(11) 


Equation  (ll)  has  been  applied  as  a  correction  to  low-speed  test  data 
for  wings  of  less  than  ±12°  sweep  and  is  compared  with  experiment  in 
figures  5  to  11;  a±  is  evaluated  from  the  low-speed  data  by  method  1. 
This  application  is  denoted  method  3. 
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Method  4.-  The  three-dimensional  Prandtl  transformation  can  he 
extended  to  thin,  swept  wings  by  means  of  Weissinger*s  approximate 
lifting-surface  theory  of  reference  4,  as  explained  in  reference  3* 
Lift— curve  slopes  calculated  hy  Weissinger,s  method  can  he  obtained 
from  charts  in  reference  3  for  a  wide  range  of  sweep  angles,  aspect 
ratios,  and  taper  ratios.  Because  the  Weissinger  method  calculations 
were  obtained  for  a  section  lift-curve  slope  of  2*  (the  thin-eirfoil— 
theory  value),  Weissingerfs  method  was  not  applied  as  a  correction  for 
compressibility  to  low— speed  data.  The  calculated  values  are  compared 
directly  with  experiment  in  figures  5  to  19  and  are  denoted  method  4. 


RESULTS  AND  DISCUSSION 
Comparison  of  Experiment  and  Calculation 


Four  methods  have  been  used  to  predict  the  variation  of  lift— curve 
slope  with  Mach  number  for  finite— span  wings.  Figures  5  to  19  show 
comparisons  between  these  methods  and  experimental  data  obtained  from 
wind-tunnel  tests  of  wing  models.  Table  I  contains  supplementary 
information  regarding  the  models  and  test  conditions  for  the  experimental 
data;  this  information  was  obtained  from  references  8,  9,  12,  13,  l4, 
and  unpublished  data.  Table  I  indicates  that,  for  models  4  to  18,  jet- 
boundary  corrections  were  either  applied  or  were  negligible.  All 
corrections  were  applied  in  a  similar  manner  to  account  for  the  effects 
of  boundary-induced  velocity  and  blockage.  Failure  to  apply  corrections 
to  the  data  of  models  1,  2,  and  3  resulted  in  a  slightly  excessive 
increase  of  lift— curve  slope  with  Mach  number  at  high  subsonic  Mach 
numbers.  Application  of  the  correction  would,  in  general,  improve 
slightly  the  agreement  between  theory  and  experiment. 

Examination  of  figures  5  to  19  shows  that  all  wing  tests  and  all 
four  methods  of  calculation  yield  an  increase  of  lift  with  Mach  number. 
Theory  and  experiment  are,  therefore,  in  qualitative  agreement. 
Quantitatively,  however,  the  agreement  between  calculation  and  experiment 
is  not  quite  consistent.  Such  inconsistencies  (figs.  8  and  10,  for 
example)  are  to  be  expected,  however,  because  all  the  calculations  are 
based  on  potential  flow.  Methods  1,  2,  and  3  are  applied  as  corrections 
to  low— speed  test  data  which  account  for  the  low— speed  boundary— layer 
effects.  If  the  Reynolds  number  changes  which  usually  accompany  Mach 
number  variation  produce  variations  in  the  boundary  layer  of  the  airfoil, 
an  additional  variation  of  lift— curve  slope  with  Mach  number,  which  is 
not  predictable  by  the  methods  discussed  herein,  will  occur.  The  marked 
difference  between  experiment  and  results  calculated  by  method  4  in 
figure  6  results  from  the  wing  model  having  a  full— span  elevator  with 
an  overhang  and  an  open  gap.  When  the  variation  of  lift— curve  slope 
with  Mach  number  was  calculated  by  methods  1,  2,  and  3,  the  gap  was 
considered  in  terms  of  its  effect  on  the  low— speed  section  lift— curve 
slope.  No  consideration  of  the  gap  was  made  when  method  4  was  applied. 
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In  general,  of  the  four  methods,  methods  1  and  2  compare  most 
favorably  with  experiment.  Of  these  two  methods,  method  1  which  includes 
a  consideration  of  section  thickness  appears  somewhat  better  than  method  2, 
particularly  for  small  sweep  angles.  Method  2  gives  less  increase  of 
lift-curve  slope  with  Mach  number  than  method  1.  Method  3  indicates 
even  less  increase  of  lift-curve  slope  with  Mach  number  than  method  2 
and  seems  to  be  in  poorer  agreement  with  experiment.  Method  4  indicates 
an  increase  of  lift— curve  slope  with  Mach  number  less  than  method  2 
but  similar  to  method  3-  Coincidence  of  the  calculations  of  method  4 
and  experiment  at  low  Mach  numbers  was  not  obtained  as  was  the  case  with 
methods  1,  2,  and  3  because  the  calculations  were  a.1 1  obtained  for  a 
low— speed  section  lift-curve  slope  of  2 it. 

The  data  of  figures  5  to  19,  either  experimental  or  calculated, 
indicate  that  the  variation  of  lift— curve  slope  with  Mach  number  decreases 
as  the  sweep  increases.  Furthermore,  a  comparison  of  methods  1  and  2 
indicates  that  the  effect  of  thickness  upon  the  variation  of  lift-curve 
slope  with  Mach  number  decreases  with  sweep,  either  positive  or  negative. 


Practical  Considerations  for  Method  1 

Because  the  agreement  between  experiment  and  method  1  seemed  good, 
a  practical  calculation  procedure,  based  on  the  assumption  that  a  low- 
speed  lift-curve  slope  is  available  for  the  finite-span  wing  under  con¬ 
sideration,  is  briefly  set  forth  for  this  method.  The  constant  apE1^  in 
equation  (10)  can  be  solved  for  as  follows: 


CLa  ^e^K 

aiFA  =  - =  K  (12) 

A  -  CT  2±i£ 

ha  it 

where  is  the  value  measured  at  some  low  Mach  number.  The  quantity  AEe 

is  evaluated  from  figure  4;  whereas  Pg  is  evaluated  from  equation  (l) 

with  being  used  rather  than  \i;  that  is,  only  the  component  of  the 

Mach  number  perpendicular  to  the  wing  quarter-chord  line  is  considered. 

The  value  of  K  from  equation  (12)  is  an  effective  section  lift— curve 
slope  for  zero  Mach  number  appropriate  to  the  wing  for  which  the  low- 
speed  data  were  obtained. 

The  finite-span  lift-curve  slope  and  its  variation  with  Mach  number 
for  the  wing  for  which  K  was  evaluated  can  now  be  written 


* 


AEe 


PK 

JL  57.3 
Pk  n 


(13) 
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If  no  low— speed  lift— curve  slope  is  available  from  which  to  evaluate  K, 

Weissinger*s  method,  which  gives  CT  at  M  =  0  (charts  of  reference  3), 

a  , 

can  he  used.  If  the  wing  is  effectively  unswept  (sweep  angles  of  quarter- 
chord  line  between  ±12° ),  K  approaches  a^,  which  can  he  determined 

from  wing-section  tests,  many  of  which  are  presented  in  reference  15. 

The  evaluation  of  (3g  requires  that  the  value  of  X  he  known  for 
the  particular  airfoil  section  under  consideration.  This  parameter  can, 
of  course,  he  calculated  as  explained  in  reference  7*  however,  such  a 
calculation  is  often  unnecessary  because  X  has  already  been  calculated, 
as  a  step  in  the  pressure-distribution  calculation,  for  a  large  number 
of  airfoils.  Figure  20  Bhows  the  variation  of  X  with  thickness  for  a 
number  of  airfoil  series.  The  value  of  X  in  equation  (l)  should 
correspond  only  to  the  airfoil  basic— thickness  form  (airfoil  without 
camber).  Equation  (l)  is,  therefore,  independent  of  airfoil  camber, 
and  figure  20  offers  a  very  wide  coverage  of  airfoils. 

As  explained  in  reference  15  there  are  two  slightly  differing  groups 
of  airfoils  both  designated  as  6-series.  The  older  group  shows  the  low- 
drag  range  following  a  comma  after  the  number  denoting  the  chordwise 
position  of  minimum  pressure  (model  1,  table  I,  for  example).  A  more 
recent  group  of  airfoils  show  the  low-drag  range  as  a  subscript  to  the 
number  denoting  the  chordwise  position  of  minimum  pressure  (model  6, 
table  I,  for  example).  The  curves  of  figure  20  apply  only  to  the  latter 
group  of  6— series  airfoils. 


CONCLUSIONS 


A  comparison  with  experiment  of  several  methods  of  predicting  the 
lift  of  wings  in  subsonic  compressible  flow  indicates  the  following 
conclusions: 

1.  Kaplan fs  formula  for  the  effect  of  compressibility  on  the  lift 
of  wing  sections  is  In  good  agreement  with  experiment. 

2.  Semi  empirical  formulas  derived  for  finite— span  wings  agree 
with  experiment  better  than  previously  derived  theoretical -methods. 

3.  A  slightly  better  agreement  at  small  sweep  angles  results  from 
a  consideration  of  thickness  in  the  semi  empirical  formulas. 
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4.  Both,  experiments  and.  calculations  indicated  a  decrease  in  the 
variation  of  lift  with  Mach  number  for  increasing  sweep. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va. ,  August  6,  19^8 
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A/gure  /  -  Vor/af/on  of  two-d/mens/or?or/  I /ff -curve 
slope  w/th  Mac/?  number.  Mode/  /.  A  =  0°. 


A/gure  2.-  Ver/at/on  of  two-d/mensio/ial  I  iff  curve 
s/ope  w/th  Mach  number ;  Mode/  2.  A=0°. 
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7/pare  3.  -  Var/at/or?  of  two -at '/eter/s, 10 nor/  lift- curve 
c/ope  w/th  Mac/?  r?urr?her  for  three  onp/es 
of  sweep .  Mode/  3 . 
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0  ./  Z  .3  .4  ,S  .6  7  .6 
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F/gure  Sr  Var/af/0/7  of  //ft-  curve  slope  w/th  Mach 
number.  A=0°.  Mode/  4. 


MCA  TW  Wo.  1739 


21 


M 

figure  8.  -  Var/af/on  of  hft-cune  s/ope  with  Mach 
number.  A  =  32°.  A/lode/  7. 


C,  10 

L/r 


1  Lowe/- limit  test  data 
fxper/me/?t _ _ _ 


Method  / 


Method  2 


Method  3 


-  Method  4 


0  ./  2  .3  .4  .5  .6  7 

M 

F/gure  d.~  Vor/ot/on  of  hft-cun/e  slope  with 
number,  A  =  32°.  Mode/  8. 
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H/pure  Hr  Var/af/or?  of  f/ff  curve  Hope  w/th  Mac/? 
number.  A  =  H5°.  /Woc/e/  /0 . 


F/pure  /3 Varya f/on  of  f/ff  carve  dope  a  iff  f/facf? 
number.  A  =  30,0°.  Mode/  //. 
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figure  /Sr  Vc/r/af/or?  0/  f/ff -curve  c/ope  w/fh 
Mac/?  m//7?/ee.  A.  =  -45,0\  Sloa/e/  /4. 
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F/pure  /7r  Van  at/ 0/7  of  t/ff -curve  a  top  a  with  Mach 
number.  A  =47.8°.  At  octet  /6 . 
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Airfoil  thickness,  percent  chord 

ftp  tv  re  20.  -  Variation  with  o/rto//  th/c knees  of  the 
a vrfot/- thickness  parameter  for  severer/  senes 
of  A/ AC  A  a /r foils . 
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